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deposition,  were  the  sane  for  all  species.  The  experimental  data 
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PREFACE 


In  the  Fall  of  1980,  we  had  high  hopes  for  continuing  our 
long-term  study  of  energy  deposition  due  to  ion  bombardment  in  metals 
and  semiconductors.  However,  due  to  shifts  in  ONR  personnel  and  program 
priorities,  we  reformulated  our  research  goals  to  study  the  stopping  of 
ions  under  "exotic"  conditions.  We  felt  that  since  high-intensity 
pulsed  irradiation  is  technologically  close  at  hand,  this  study  appeared 
timely  and  worthwhile.  Nevertheless,  in  the  decision  process  this 
direction  was  considered  to  be  too  far  afield  from  the  program  objectives. 
Very  kindly,  ONR  did  grant  us  funding  for  what  was  to  be  an  extension  to 
ONR  N00014-76-C-0482,  in  order  to  tidy  things  up  and  write  a  final  report. 
Due  to  further  delays,  however,  it  became  impossible  to  extend  the  old 
contract  and  a  new  one  had  to  be  issued.  This  meant  that  a  final  report 
now  had  to  be  written  in  the  late  Fall  of  1931  covering  the  period 
7  January  1976  to  31  October  1981. 

This  rather  complicated  series  of  events  left  the  new  contract 
N00014-81-C-0632  with  its  main  purpose  accomplished  and  reported  upon. 
However,  not  anticipating  all  these  events  early  enough,  we  had  started 
a  research  effort  on  the  simulation  of  radiation  damage  in  semiconductors 
quite  complementary  to  our  previous  work  in  metals.  By  scraping  together 
various  resources,  we  have  been  able  to  complete  the  studies  on  single 
crystal  Si  and  give  these  results  here  as  the  main  body  of  this  final 
report.  Readers  interested  in  our  previous  work  are  referred  to  the  final 
report  of  contract  N00014-76-C-0482. 
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SUMMARY 


-1/. 


A  combined  theoretical  and  experimental  study  of  primary 
recoil  spectra  effects  or  radiation  damage  in  silicon  is  presented. 


Calculations  determined  how  the  damage  energy  is  partitioned  into  free 


defects  and  cascajdes  by^fast  collisions.  The  theory  also  showed  that 
on  a  time  scale  sec,  a  very  weak  mass  dependence  of  the  lattice 

damage  is  to  be  expected.  Channeling  experiments  were  then  performed 


on  <111>  single  crystal  silicon  implanted  with  1.0  MeV.wfte,  0.5  MeV.l  -He, 
fl\  $  i~J  #  ^ 

and  75  keV^^H.  Energies  and  fluences  of  the  ions  were  matched  such  that 


over  the  first  0.3  the  damage  energy  deposited  and  the  rate  of 

\ 

energy  deposition  were'  the  same  for  all  species.  The  experimental  data  were 
analyzed  assuming  that  equivalent  primary  damage  states  will  evolve  into 
statistically  equivalent  final  damage  states  at  high  fluences.  They 
confirm  that  the  final  damage  is  essentially  independent  of  the  mass  of 


the  bombarding  ion. 


A 


tu<  CYX> 


1 .  INTRODUCTION 


A  primary  concern  in  the  simulation  of  neutron  damage  by  ion 

beams  and  in  the  study  of  ion  beam  induced  damage  in  general  revolves 

around  how  the  initiating  (incident)  particle  interacts  with  the  target 

to  produce  recoil  atoms  and  atomic  displacements.  Some  of  the  key 

questions  relate  to  understanding  how  the  energy  is  partitioned  in  the 

-14  -12 

host,  both  on  a  fast  collision  time  scale  of  10  to  10  sec  and  on 
the  longer  time  scale  associated  with  cascade  collapse.  One  needs  to 
determine  whether  the  amount  of  damage  created  by  particles  of  different 
mass  and  energy  depends  on  the  total  energy  deposited  into  displacement 
processes.  Another  fundamental  question  pertains  to  the  spatial  distribu¬ 
tion  or  degree  of  localization  of  the  damage  produced  by  recoiling  atoms 
of  different  energies. 


Considerable  interest  in  the  above  problems  prompted  a  combined 
theoretical  and  experimental  study  of  primary  recoil  spectra  effects  at 
the  HEIBS  (High  Energy  Ion  Bombardment  Simulation)  facility  at  the  University 
of  Pittsburgh.  The  theoretical  effort  was  aimed  at  calculating  the 
relative  energy  deposited  in  free  defects  and  subcascade  regions  for 
energetic  recoils  (E  i  1  keV)  before  annealing. 


There  is  considerable  experimental  evidence  and  some  previous 

theoretical  calculations  which  suggest  that  the  energy  in  a  displacement 

cascade  does  not  increase  indefinitely  with  increasing  primary  recoil 

atom  energy  in  materials  where  cascades  can  occur.  Rather,  above  some 

energy,  E  ,  cascades  split  into  well-defined  separate  subcascade  regions. 

(1)U 

Merkle  has  summarized  much  of  the  experimental  evidence  for  subcascade 


*The  HEIBS  effort  is  jointly  comprised  of  Westinghouso  R&D  Center  and 
University  of  Pittsburgh  personnel. 
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formation,  and  their  existence  has  also  been  suggested  by  computer 

simulations  of  atomic  displacement  cascades  using  the  binary-collision 
(2  3) 

approximation.  ’  In  this  work,  the  subcascade  production 
probability  was  calculated  on  the  basis  of  the  LSS-LNS  theory  of  ion 
scattering  and  energy  loss^*  and  clearly  shows  that  the  primary 
damage  state  is  remarkably  independent  of  the  ion  species  and  its  PKA 
spectrum. 

Since  it  is  almost  impossible  experimentally  to  determine 
the  partitioning  of  energy  into  subcascades  and  "free"  defects  on  the 
fast  collision  time  scale,  experiments  were  designed  cn  the  premise  that 
statistically  equivalent  primary  damage  states  (calculated)  will  evolve 
into  statistically  equivalent  final  damage  states  at  high  finances. 

This  means  that  techniques  used  to  "measure"  the  resulting  damage,  such 
as  Rutherford  backscattering/channeling  experiments,  should  not  be  able 
to  distinguish  between  damage  resulting  from  ions  of  different  mass  and 
energy  (and  therefore  PKA  spectra)  if  the  total  damage  energy  deposited 
into  the  target  and  the  rate  of  energy  deposition  are  the  same  for  all 
ions.  This  presupposes  that  the  ions  are  sufficiently  energetic  to 
launch  recoils  which  can  produce  cascades.  To  determine  whether 
differences  in  the  primary  recoil  energy  spectra  result  in  measurably 
different  damage  states,  single  crystal  silicon  targets  were  bombarded 
with  ions  of  different  mass  and  energy  and  the  resulting  damage  states 
compared.  This  comparison  was  accomplished  by  Rutherford  backscatterinb/ 
channeling  studies  of  the  damaged  crystals. 
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2.  PKA  SPECTRA  AND  CASCADES 

Primary  knock-on  atoms  (PKAs)  are  atoms  of  the  target  substance 
which  directly  receive  energy  in  a  collision  with  the  initial  radiation 
(neutrons,  ions,  or  electrons).  The  number  of  PKAs  as  a  function  of  their 
energy  constitutes  the  PKA  speetrumi  it  is  directly  determined  by  the 
collision  kinematics  and  the  differential  cross  section.  For  example,  the 
maximum  PKA  energy  is  determined  by  the  masses  of  target  and  incident 
particle. 

Electron  beam  irradiation  produces  only  low-energy  PKAs  because 
cf  the  great  mass  difference  of  electrons  from  target  atoms.  Much  more 
energetic  PKAs  are  produced  by  fast  neutrons  (E  £  1  MeV),  but  because  the 
neutrons  are  uncharged  they  do  not  produce  very  many  low-energy  PKAs. 
Although  the  ion  PKA  spectrum  can  be  controlled  to  some  extent  by  selection 
of  ion  species  and  energy,  it  is  clear  that  one  cannot  exactly  match  a 
realistic  neutron  PKA  spectrum  with  fast  ions  since  the  ions  will  always 
produce  many  more  low-energy  PKAs  than  will  neutrons. 

The  dissimilarity  of  PKA  spectra  naturally  leads  to  concerns 
about  the  suitability  of  ion  beams  to  simulate  neutron  damage  in  materials. 
At  first  sight,  it  appears  that  for  ions  a  harder  PKA  spectrum  will  cause 
more  energy  to  be  deposited  in  cascades  relative  to  free  defects.  Careful 
examination  of  the  subcascade  production  process  reverses  this  impression. 

PKA  events  have  been  divided  into  three  categories  depending  on 
the  recoil  energy: 

1.  low-energy  PKAs  (below  1-2  keV)  which  produce 
defect  pairs  or  only  small  clusters; 

2.  intermediate-energy  PKAs  (2  -  12  keV)  which  produce 
single  cascades;  and 

3.  high-energy  PKAs  (E  =  12  keV)  which  produce 
multiple  subcascades  and  free  defects. 
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The  boundary  energies  are  not  sharply  defined,  and  calculations 

have  been  performed  to  explore  the  range  of  uncertainty  in  the  boundary 

energies. ^  The  probability  calculations  strongly  support  the  general 

division  into  three  energy  categories;  for  example,  our  calculations  show 

that  an  18  keV  PKA  in  Si  has  a  probability  greater  than  90%  to  launch  a 

subcascade  before  dropping  below  'v 12  keV  to  form  its  own  subcascade.  For 

this  reason,  that  part  of  the  PKA  spectrum  which  lies  above  Eu  'v  12  keV 

has  no  direct  physical  significance  or  manifestation;  PKAs  in  this  range 

are  instead  converted  into  lower-energy  subcascades  and  free  defects. 

-14 

This  conversion  occurs  on  the  fast-collision  time-scale  ('v  10  sec) 
and  is  complete  before  any  thermally  activated  atomic  motions.  Because  of 
this  conversion,  the  hard  neutron  spectrum  produces  free  defects  and 
subcascades  very  much  like  the  damage  resulting  from  ion  irradiation. 

Because  the  process  discussed  here  is  relatively  energetic 

scattering  (E  £  12  keV),  it  can  be  modeled  with  reasonable  accuracy 

by  the  screened  Coulomb  potential  introduced  by  Lindhard,  Neilsen  and 
(4) 

Scharff.  Calculations  of  the  probability  to  launch  subcascades  and 
the  demonstration  that  an  energy  Eu  exists  such  that  subcascades  are 
almost  certain  to  be  launched  for  PKA  energies  above  Ey  have  been 
described  in  detail  previously .  ^  Similarly,  this  refernece^  also 
discusses  the  specifics  of  the  sequel  calculations  which  follow  the 
collision  sequence  in  which  a  PKA  having  energy  greater  than  Eu  is 
allowed  to  transfer  energy  to  electrons,  low-energy  recoils,  and 
subcascade-forming  recoils.  What  will  be  presented  here  are  the  results 
of  such  calculations  as  performed  for  self-ion  bombardment  of  silicon. 
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3.  EXPERIMENTAL  BACKGROUND 


3.1  Selection  of  Target  Materials 

Single  crystal  specimens  have  been  selected  to  permit 

application  of  the  Rutherford  backscattering/channeling  technique.  The 

Si  specimens  were  cut  from  crystals  with  a  <111>  growth  axis.  Small 

wafer  specimens  of  approximately  1.5  cm  on  an  edge  and  250  pm  thick  were 

cleaved  along  {111}  directions  from  larger,  250  pm  thick,  wafers. 

15  3 

Impurity  levels  in  this  material  are  ^  10  /cm  . 

3.2  Choice  of  Bombarding  Ions  and  Implant  Conditions 

ince  the  experiments  utilized  ions  of  different  mass,  both 

the  energy  of  the  bombarding  ions  and  their  flux  had  to  be  selected  such 

that  the  total  damage  energy  deposited  into  the  target  and  the  rate  of 

energy  deposition  was  the  same  for  all  ions.  In  addition,  chemical 

effects  from  the  implanted  ions  and  from  point  defect  concentration 

gradients  induced  by  proximity  of  a  free  surface  had  to  be  minimized. 

Hence,  in  this  study,  the  inert  gas  ions  ^Ne  and  ^He  as  well  as 

20  + 

were  employed.  The  1  MeV  Ne  ions  were  selected  as  the  reference  ion 
and  energy  for  establishing  the  damage  profile,  SD(x)  in  the  Si  crystals. 
Fluences  and  energies  for  the  other  two  species  were  then  determined 
to  provide  as  close  a  match  as  possible  to  the  reference  Ne  damage 
profile.  There  are  difficulties  with  the  present  range-energy  theories 
but  we  believe  that  S^(x),  energy  oer  unit  depth  deposited  into  atomic 
displacements,  may  be  calculated  reliably  for  the  region  starting  at  the 
surface  and  going  into  the  sample  to  a  depth  corresonding  to  approximately 
three-quarters  of  the  way  to  the  peak  of  the  displacement  damage.  Figure  1 
shows  curves  of  S^(x)  calculated  from  a  modified  EDEP-1  code  of  Manning 
and  Mueller^  for  75  keV  ^H,  0.5  MeV  ^He,  and  1  MeV  ^Ne,  which  yield 
equal  energy  deposition  near  the  surface  of  the  Si  target  (0  to  0.3  pm). 


x) ,  MeV/p  m/atom 


The  Sr^x)  curves  are  normalized  at  the  front  surface  and  the  fluences 
adjusted  to  agree  with  20Ne  damage. 


Table  1  presents  the  choice  of  bombarding  ions,  energies,  and 
relative  fluences  selected  to  produce  approximately  matching  damage  profiles 
over  distances  of  ^0.3  ym.  Also  included  are  the  relevant  range  parameters 


(Rp  and  A R  )  for  each  ion,  and  the  maximum  transfer  energy,  T^ 


Table  1  -  Experimental  Parameters 


Energy 

Ion 

(MeV) 

20Ne 

1.0 

*He 

0.5 

h 

0.075 

(pm) 

(ym) 

1.52 

0.168 

1.95 

0.122 

0.71 

0.063 

T 

max 

(keV) 

1 

972 

38 

219 

99 

9.98 

4.  EXPERIMENTAL  PROCEDURES 

The  final  silicon  single  crystal  surfaces  utilized  for 
implantation  were  perpendicular  (within  1  to  2°)  of  the  [111]  growth 
direction  and  were  free  from  mechanical  damage.  Specimens  were  implanted 
5  to  7°  off  the  [111]  axis  normal  to  the  surface  to  minimize  channeling 
effects  and,  during  the  bombardment,  one-half  of  each  was  masked  to 
provide  a  nonimplanted  reference  crystal  for  the  channeling  studies. 

The  channeling  experiments  were  performed  on  a  4.8  meter  beam 
line  of  the  2  MV  Van  de  Graaff  accelerator  using  1.5  MeV  ^He+  ions.  The 
^He  beam  was  collimated  to  give  a  full  angular  divergence  of  0.03°.  For 
the  aligned  spectra,  the  He  beam  was  oriented  normal  to  the  crystal  surface 
(within  2°)  and  all  spectra  were  collected  at  a  backscattered  angle  of 
168°.  The  aligned  spectra  were  obtained  in  a  <111>  axial  direction  from 
both  the  implanted  and  nonimplanted  half  of  each  crystal  and  the  random 
spectrum  was  obtained  from  the  nonimplanted  half.  The  random  spectrum 
was  measured  by  setting  the  symmetry  axis  8°  from  the  incoming  beam 
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direction  and  rotating  the  sample  continuously.  Each  spectrum  was 

obtained  under  a  constant  total  He  fluence  as  determined  by  the 

backscattered  yield  from  a  rotating  Au  foil,  which  sampled  the  beam 

approximately  7%  of  the  time.  The  energy  scale  of  the  backscattered 

spectra  was  converted  to  a  depth  scale  by  employing  the  electronic 

4  + 

stopping  powers  of  He  in  Si,  with  the  assumption  that  the  energy  loss 

of  the  channeled  ions  is  the  same  as  that  for  ions  impinging  along  a 
(9) 

random  direction. 

5.  RESULTS  AND  DISCUSSION 

5.1  Subcascade  Spectra 

In  the  calculations,  it  was  considered  that  a  cascade  was 
produced  by  a  recoil  of  energy  greater  than  ^0.6  -  2  keV.  As  E^ 
is  varied  through  this  range,  the  upper  limit  Ey  (above  which  multiple 
cascades  are  produced)  will  also  vary.  However,  the  general  features  of 
the  energy  partitioning  remain  constant. 

Figure  2  presents  the  calculated  direct  and  reduced  energy 
spectra  for  various  PKA  energies  in  self-bombardec  silicon.  The  direct 
spectrum  shows  the  energy  deposited  in  secondary  events  by  PKAs  of 
specified  energy.  Values  of  E^  and  Eu  were  1.6  keV  and  12  ksV, 
respectively.  As  the  PKA  energy  rises,  a  larger  fraction  of  tha  energy  is 
delivered  to  energetic  (E  >  Ey)  secondaries  (shaded  box) .  The  reduced 
spectrum  is  the  result  of  distributing  the  energy  of  the  shaded  box  into 
energy  of  free  defects,  cascades,  and  electrons.  A  most  important  finding 
is  that  the  partitioning  of  the  energy  into  tree  defects  and  cascades,  as 
indicated  by  the  reduced  spectra,  is  almost  the  same  for  all  recoil 
energies  up  to  the  maximum  energy  transfer  permitted  by  kinematics.  A 
schematic  interpretation  of  v'-at  these  results  mean  in  terms  of  defect 
and  cascade  formation  is  shown  in  Figure  3.  Since  the  ratio  of  reduced 
spectra  histogram  heights  in  Figure  2  is  almost  independent  of  PKA  energy, 
it  is  relatively  insensitive  to  the  PKA  spectrum. 
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Figure  2.  Calculations  of  the  direct  and  reduced  energy  spectra  for 

various  PKA  energies.  The  direct  spectrum  shows  the  energy 
deposited  in  secondary  events  by  PKAs  of  specified  energy. 

'll  is  1.6  keV  and  Eu  is  12  keV  in  this  figure.  As  the  PKA 
energy  rises,  a  larger  fraction  of  the  energy  is  delivered  to 
energetic  (E  >  Eu)  secondaries  (shaded  box).  The  reduced 
spectrum  results  from  distributing  the  energy  of  the  shaded  box 
into  energy  of  free  defects,  cascades  and  electrons.  In  the 
reduced  spectrum,  the  cascade  energy  is  a  constant  multiple 
of  the  free  defect  energy,  independent  of  PKA  energy. 


5.2  Backscattering  Data 

Figures  4  through  6  show  aligned  backscattered  spectra  for 

1  MeV  ^Ne.  0.5  MeV  ^He,  and  75  keV  ^H,  respectively.  Also  shown  on 

each  figure  is  a  representative  aligned  spectrum  for  a  nonimplanted 

(n-i)  reference  crystal  and  a  random  nonaligned  spectrum.  All  yields 

are  normalized  with  respect  to  the  random  spectrum.  The  peak  at  zero 

depth  for  the  aligned  spectra  is  due  to  the  direct  backscattf  '-ing  from 

the  surface  layer  of  ions  which  do  not  enter  a  channel  and  is  known  ?.s 

the  "surface  peak."  For  the  bombarding  ions,  the  surface  peak  is 

essentially  the  same  for  1.39  x  10^  Ne+/cm^  and  5.31  x  10^  He+/cm^, 
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but  is  somewhat  larger  for  1.37  x  10  H  /cm  ,  perhaps  due  to  the 
increase  in  fluence.  The  "errors"  in  the  backscattered  spectra  are- 
represented  by  the  recorded  yield  variations  on  the  curves  (i.e.,  the 
small  peaks  and  valleys). 

Comparison  of  the  implanted  spectra  behind  the  surface  peak 
reveals  very  similar  behavior  up  to  depths  of  ''<400  nm,  or  slightly  beyond 
the  region  over  which  the  deposited  damage  energy  was  matched  (Figure  1). 
(Note  that  the  large  peak  at  ''<600  nm  in  Figure  6  corresponds  to  the  peak 
of  the  SpCx)  curve  at  640  nm.)  The  spectra  shown  in  Figures  4  and  6  can 
be  compared  more  easily  at  equivalent  damage  levels  by  an  approach 
outlined  by  Merkle  et  al.^*^  This  analysis  is  frequently  applied 
in  channeling  experiments  and  relates  the  rate  of  dechanneling  to  the 
concentration  of  defects  of  type  (j)  and  density  (n^)  with  a  cross 
section  for  dechanneling  par  defect  (cjj).  If  one  normalizes  the  aligned 
spectrum  by  the  random  spectrum,  the  normalized  yield  (x)  represents  the 
dechanneled  fraction  of  ions  up  to  the  depth  (z),  and  (1-x)  represents  the 
remaining  channeled  fraction.  Therefore,  at  any  depth  z: 
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Figure  A.  Comparison  of  normalized  yield  spectra  (<111>  axis)  for 
silicon  implanted  with  1.0  MeV  Ne  to  a  fluence  of 
1.39  x  10^  ions/cm  and  nonimplanted  silicon  (curve 
marked  (n-i)).  Implants  were  made  at  approximately  5° 
off-axis  and  all  spectra  were  obtained  with  a  1.5  MeV 
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Figure  5.  Comparison  of  normalized  yield  spectra  (<lll>  axis)  for 

silicon  implanted  with  0.5  Hey  He+  to  a  fluence  of  5.31  x  10 
ions/cm2  and  nonimplanted  silicon  (curve  marked  (n-i)). 
Implants  were  made  at  approximately  5°  off-axis  and  all 
spectra  were  obtained  with  a  1.5  MeV  4He  beam. 
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Figure  6.  Comparison  of  normalized  yield  spectra  (<111>  axis)  for 

silicon  implanted  with  75  keV  H+  to  a  fluence  of  1.37  x  1017 
ions/cra2  and  nonimplantea  silicon  (curve  marked  (n-i)). 
Implants  were  .  e  at  approximately  5°  off-axis  and  all 
spectra  were  ouuained  with  a  1.5  MeV  4He  beam. 


where  the  subscript  (i)  refers  to  the  aligned  spectrum  from  the  implanted 
crystal.  This  equation  simply  states  that  the  rate  of  dechanneling 
at  any  depth  is  given  by  the  fraction  of  ions  remaining  in  the  channels 
times  the  probability  that  a  single  ion  will  be  dechanneled  by  all  the 
defects  of  type  (j)  as  it  progresses  a  unit  distance  along  the  channel. 

If  the  intrinsic  dechanneling  as  measu.  ^d  in  the  nonimplanted  reference 
crystal  is  simply  additive  to  that  produced  by  tne  defects,  then  equation  1 


can  be  written  as 

% 


where  the  subscripts  (i)  and  (n-i)  refer  to  the  spectra  from  the  implanted 
and  nonimplanted  crystals,  respectively.  This  expression  strictly  holds 
when  the  defects  introduce  dechanneling  of  the  ions  into  random 
trajectories  without  any  direct  backscattering  of  ions  in  the  channels. 
Therefore,  according  to  equation  2,  the  density  of  defects  and  thus  the 
damage  level  as  measured  by  a  channeling  experiment  is  related  to  the 


slope  of  the  aligned  spectrum. 

k 


The  quantity 
which  first  determines 


I 

j=l 


o^n^  was  calculated  using  a  computer  program 
the  slopes  of  the  channeling  spectra.  The 


depth  range  was  selected  to  avoid  the  proximity  of  the  small  surface 


peak  and  to  extend  only  to  the  limit  over  which  the  SD(x)  curves  were 
matched. 


Figure  7  shows  a  plot  of  E  a.n.  as  a  function  of  depth  (as 

J  3  15 

determined  by  equation  2)  for  the  damage  level  equivalent  to  1.39  x  10 

Ne  /cm  .  The  data  are  also  tabulated  in  Table  2.  It  is  clear  that  the 

extent  of  dechanneling  for  the  damaged  crystals  bombarded  with  1.0  MeV 
20  4 

Ne  and  0.5  MeV  He  superimpose,  which  suggests  that  equivalent  final 
damage  states  were  produced  for  these  two  ions  This  further  suggests 
that,  provided  the  deposited  damage  energy  and  the  rate  of  energy  deposition 
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are  the  same,  the  damage  state  produced  is  essentially  independent  of  the 

mass  and  energy  of  the  bombarding  ion  in  the  silicon  substrate.  Such  a 

conclusion  is  in  excellent  agreement  with  data  recently  reported  on 

(12) 

ion-bombarded  molybdenum. 


Table  2  - 

Dechanneling  Rates  as  a 

Function  of  Depth 

Implanted  Ion 

Depth  (ym) 

Z  o.n.  C1, 

1  l 

20 

1.0  MeV  Ne 

0.11 

0.14 

0.22 

0.18 

0.31 

0.19 

0.41 

0.19 

0.5  MeV  4He 

0.10 

0.14 

0.20 

0.15 

0.30 

0.16 

0.40 

0.15 

75  keV  1H 

0.13 

0.23 

0.20 

0.33 

0.30 

0.48 

0.39 

0.62 

The  higher  slope  of  the  75  keV  ^H 

curve  corresponding  to 

a  greater  dechanneling  rate,  is  due  to  the  influence  of  the  large  peak 


observed  at  0.6  ym  in  Figure  6.  The  implanted  hydrogen,  with  Rp  =  0.71  ym, 
is  introducing  lattice  strain  which  yields  an  additional  contribution 


to  the  fraction  of  dechanneled  ions  in  the  near  surface  region.  This 

does  not  necessarily  indicate  that  the  damage  state  produced  by  the  75  keV 

^H  as  it  enters  the  crystal  differs  from  the  states  produced  by  1.0  MeV 
20  4 

Ne  and  0.5  MeV  He.  Further  work,  including  electron  microscopy,  is 
required  to  elucidate  this  question.  One  might  expect  that  the  resultant 
damage  state  for  ^H  would  differ  from  those  produced  by  ^Ne  or  Sie 
because  the  recoil  spectrum  for  75  keV  in  silicon  would  primarily 
yield  free  defects.  However,  since  the  maximum  transfer  energy  is  still 
'vlO  keV  in  silicon,  which  is  considerably  larger  than  E^,  a  significant 
amount  of  cascade  formation  could  still  occur. 
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An  interesting  question  which  is  still  to  be  answered  concerns 

the  nature  of  the  damage  produced  by  the  ion  bombardment.  It  is  hoped 

that  transmission  electron  microscopy  (TEM)  will  be  able  to  determine 

whether  amorphous  regions  have  been  produced,  or  if  the  dechanneling  is 

(13) 

primarily  due  to  loop  formation,  for  example.  Ligeon  has  suggested 

that  30  keV  H  implanted  into  silicon  does  not  produce  amorphicity, 

even  at  the  peak  damage  region,  but  that  hydrogen  bubbles  may  be 

contributing  to  the  observed  dechanneling.  This  is  difficult  to 

rationalize  because  of  the  weak  dechanreling  effect  normally  caused  by 

(13) 

bubbles,  and  so  Ligeon  further  suggests  that  the  bubbles  are  active 
as  generalized  stacking  faults.  TEM  characterization  of  the  implanted 
region  nay  shed,  light  on  this  question,  as  it  has  in  previous  work 
on  molybdenum. 


6.  CONCLUSIONS 

•  Theoretical  calculations  have  shown  that,  on  the  fast-collision 
time  scale,  the  ratio  of  energy  deposited  into  free  defects  and 
subcascades  is  independent  of  PKA  spectrum  and  thus  of  incident  ion 

energy  and  mass  (provided  the  damage  energy  deposited  and  the  rate  of 
energy  deposition  remain  the  same). 

•  Channeling  experiments  on  <111>  single  crystal  silicon 

20  +  4  x  14. 

implanted  wi  :h  1.0  MeV  Ne  ,  0.5  MeV  He',  and  75  keV  H  ions  have 

shown  that  equivalent  final  damage  states  were  produced  for  all  three 

ions . 

•  These  two  results  further  support  the  premise  that  equivalent 
primary  damage  states  will  evolve  into  statistically  equivalent  final 
damage  states  at  high  fluencos. 
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